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Abstract 

Exitation of atomic levels due to interaction with electromagnetic waves has been the subject of numerous 
works, both experimental and theoretical. This topic became of interest in accelerator physics in relation to 
high efficiency charge exchange injection into rings for high beam power applications. Taking equations of 
resonant atom-wave interaction equations as a basis, this paper shows that there exist some interesting 
phenomena which lead to the existence of trapped electomagnetic waves (photon traps) in a medium that 
consists of atoms with transition frequencies in proximity to the wave frequency. These traps may exist in 
random and periodic lattices, and may have very low loss rate. The atomic medium can serve as an 
excellent wavegiude or tool to form and transmit electromagnetic waves for applications to accelerators and 
to electromagnetic devices in general, where high pressure gas use is acceptable. In addition, such traps in 
gases may accumulate substantial energy for a long period of time, leading to the possibility of creating 
objects similar (or equivalent) to ball lightning. 

I. INTRODUCTION 

Periodic lattices of atoms, placed in a monochromatic electromagnetic wave with its 
frequency close to that of some transition frequency between atomic levels, produce either 
oscillator}' solutions for the wave function in the space domain, corresponding to 
propagating waves, or exponentially growing/ decaying solutions in the limit of zero electric 
field. When the fields are strong, nonlinearity changes the nature of the interaction, giving 
rise to finite in space (or trapped) solutions for the coupled atom-wave equations. The 
equation for trapped waves in the Rotating Wave Approximation (RWA) (see, e.g., [1]) for a 
ID lattice of simple hydrogen-like atoms is given in the next section. Section III introduces 
approximate solutions for these equations. Appendix A provides trapped wave numerical 
solutions for ID random atom lattices. It is shown that the localization phenomenon, 
discovered by Anderson in connection with random lattices in metallic alloys (see [2]) may 
play a significant role in creation of trapped modes as well as in the nonlinearity of the atom- 
wave interaction. Section IV presents a generalization of the results for 2D and 3D cases for 
an averaged equation in the case of high atom density. Section V gives practical estimations 
for the achievable electric field of the trapped wave in the case of real media. The conclusion 
summarizes the results. 

II. ID SEMICLASSICAL EQUATIONS FOR RESONANT ATOM- WAVE 

INTERACTION 

To obtain the solution for trapped electromagnetic fields, the easiest approach is to use semi 
classical equations for two level atoms, interacting with a monochromatic classical 
electromagnetic wave. The atoms may have more than two levels, but we consider here only 
the resonant case in which some two level transition frequency C0 is close to that of the 
wave frequency CO. In this section, we deal also with a periodic ID lattice, consisting of 
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straight sections, where the electromagnetic wave propagates freely, and of infinitely short 
2D planes, consisting of atoms with area density N arranged perpendicular to the wave 
propagation direction. The equation for slow varying* probability amplitudes C 1 (lower 
level) and C 2 (upper level) is (see [1]): 

• _ i/l l2 E m 

• _ l/J, 2l rL Q A , 

2 2h 1 

where A=0>0) (h the electric field has the form E=E cosOM, 

jU l2 = jU* 2l = —J d 3 r u[ (r)ezu 2 (r) (assuming the light is polarized in the direction % parallel 

with the atomic plane), and Uj and u 2 are the normalized wave functions of the lower and the 
upper excited states, respectively. In the case of hydrogen, the lower level has primary 
quantum number n= 1 . The upper level primary quantum number is determined by the 
resonant condition, its angular momentum 1=1, and the projection of angular momentum on 
the z axis m=0. For this case jU=jU ltl is real and we omit its subscripts for simplicity. 

These equations have an interesting solution in which both amplitudes C, and C 2 oscillate at 
fixed frequencies, with no change of their absolute values: 



C l2 =\C l2 \ exp(/<Wj 2 t). (2) 



Substituting (2) to (1) yields: 



co l =co 2 + A, 

-a±J¥7q y ( 3 ) 

«*= 2 

where Q = juE / h is the Rabi frequency at exact resonance. It is clear that the atomic 

excitation is constant in time for these solutions and that the electric current, caused by 
oscillations of probability amplitudes, has a 90 degree phase shift with respect to the electric 
field. Thus, the energy of the field is neither absorbed nor emitted by the atoms. To calculate 
atomic reflection and transmission of the electromagnetic wave, we use formulas for the 
atomic plane polarization P induced by the electric field, and Maxwell's equations for the 
electric field, influenced by atomic currents, from [3]: 



a Slow varying amplitudes means that the real amplitudes are slow amplitudes, multiplied by exp(-iG))t), where 0) u 
is the corresponding level eigenfrequency. 
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P z (x,t) = 2NjuReC lf C 2f *S(x), 
d 2 d 2 _ 1 d 2 P z (x,t) ( 4 ) 

\ 2 2^ 2^ z 2 ~\ 2 ' 

ox c dY £ c at 

where N is the area density of atoms, the atomic plane is located at x=0, and all atoms have 
same frequency. The subscript/ for probability amplitudes means we have to use their "fast" 
values and real frequencies of levels (see footnote on page 2). Now we are in position to 
calculate the transmission and reflection coefficients of the wave at the atomic plane. The 
electric field atx=0 is E=EgCosOX, and the solution of (1) for "slow" probability amplitudes 
becomes: 




a 2 



(5) 

exp(ia) 2 t). 



Multiplying by factors exp(-iO) l1J2 t), their own Eigen frequencies, and substituting into (4), 
the polarization P K (x,t) becomes: 



n, x Nf£l , , Nju z E(x = 0) . 

P z (x,t) = A ^ 2 cos(0)t)S(x) = / S(x), (6) 

e a) 2 (l + (—^)) e ha) 2 (l + (—^)) 

Aco 2 Aco 2 

where <»is the electromagnetic wave frequency. One can see that the polarization oscillates 
either in phase or anti phase with electric field, depending on sign of C0 2 . The current, which 
is proportional to the polarization derivative, oscillates with a phase shift of 90 degrees 
relative to the field phase. To calculate the transmission and reflection coefficients, we take 
the electric field coming from one (positive x) side of the atomic plane as a sum of incident 
Ej and reflected E r waves. On the other side of the plane we have, for simplicity, only the 
transmitted wave E t . In mathematical notation, they are: 

E i+ = E cosk(x + ct), 

E r+ = E l sin k(x - ct) + E 2 cos k(x - ct), (7) 
E t _ = E 3 cos k(x + ct) + E 4 sin k(x + ct), 

where E t are the amplitudes of the waves, k=Q)/ c, and the signs + and - stand for positive 
and negative x coordinate. From equations (4) it follows that the electric field is continuous 
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, . i • • r , , • Nk 2 u 2 E(x = 0) 
at x— t/, and the derivative or the electric field with respect to x jumps 



e hoj 2 (l + (—-j)) 
Aa> 2 



•>2 



This yields: 



Nk 2 jU 2 E 

Aa% 
Nk 2 ju 2 E, 



(E l -E 4 )k= ^2 ' (8) 

£ /z<» 2 (l + ^— r ) 



where <W^i2 should be taken at the electric field amplitude -J is-, 2 + £4 which is the modulus 
of the electric field at x=0. 

Equations (7-8) can be rewritten in more compact complex form: 

E f _ + E b _ = E f + + E b+ , 

iNkjU 2 (E f +E h ) 
(E f+ ~E f _) = ^ W J , (9) 

2f ^« 2 (l + — ) 
Aco\ 

where the complex amplitudes of the forward and backward waves are related to real fields 
through E real = Re E f(b) exp(—ik(ct + x)) b . The subscripts / and b denote the forward and 

backward directions (forward means the direction to +infinity), and the subscripts + and - 
are related to the fields immediately to the right and to the left of the atomic plane, 
respectively. One can check that equations (9) are valid when we take the incident wave in 
(7), coming from the other direction. When one moves one step in the positive direction 
from one lattice node to another, the complex amplitudes of the forward and backward 
waves are multiplied by factors exp(— and exp(/0) , respectively, where 0is the wave 
phase advance between the two atomic planes. Combining these factors with (9), we obtain 
the transformation of the wave complex amplitudes over one lattice period: 



b The sign of factor in front of / in the exponent should be the same for forward and backward waves to satisfy 
(8) and (9). 



4 



iNkjU 2 (E fn +E hn ) 
E f, n+ i = expH0)(£ />n 2 ' ), 



2£ () tlC0 2 (\ + - r7 j) 

(10) 



A(o l 2 



iNkju 2 (E fn +E hn ) 
E h , n+ i = exp(/0)(£ 6jB + 2 ' ), 

2eji(a 1 (\ + — T ) 
A(o\ 

where we must take (0 2 ,Q at the electric field amplitude I E f _ + E b _ I . The index n stands for 

the lattice node with ascending numbers corresponding to the "right" direction. It is 
convenient to normalize the electric field such that its Rabi frequency is expressed in terms 

2uE 

of frequency offset A. The new variable for this normalization is u = . This transforms 

HA 

the equations (10) to the following form: 

_ ifc(u fn +u hn ) 
u f ,„ + i =expH0)O /n + — = 



lK{u fn +U bn ) 

u b , n+ i = expO'0)(w M ± i J ' —), 



2 Jl+ I U f" +U »- | 2 
V 2 (n) 



2jl + \ Uf " +Ub ' n I 2 



I// 2 

where K = 



£ Q hA 



To find stationary solutions, we can further simplify this map when we have u^c.c. u h in 
some straight section. One can check that the map (11) preserves this relation of the waves: 
the forward wave complex amplitude will always be the complex conjugate of the backward 
wave amplitude. Taking Re u^q, Im u^p we have: 

K ' n 

<ln + X = In C0S + iPn + I " , ) Sin 0> 

v (12) 

Pn +l = 'In Sln ^ + (Pn + ) " ) COS </>■ 

Map (12) is a typical ID symplectic map of classical mechanics. The phase space of this 
motion has stable and unstable points, closed invariant curves, chaotic trajectories, etc. Here 
we concentrate mosdy on the behavior of the fields near the phase space center point with 
coordinates (0,0). 

For small K~1 the behavior of the solution near zero is of two types: stable motion 
around the center, and unstable motion, which approximately follows a closed curve away 
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from and returning to the coordinate center (0,0). Figure 1 shows the map's (12) (with the 
sign chosen negative in the first and second equation) electric field behavior over the lattice 
nodes with initial conditions (x,p)= (0.00002, 0.00002) and parameter K=-0.5 and 0=1.0 
radian. When the parameter K increases, the center point (0,0) becomes unstable. Linearizing 
(12) around zero, and finding eigenvalues of the matrix, one gets the condition for stability: 

Icos0 + — sin^k 1. (13) 




FIG. 1. Electric field oscillation in ID lattice, consisting of resonant atom planes. Horizontal 
and vertical axes show real and imaginary part of the forward field, respectively. Parameter K 

=-0.5. 

Figure 2 shows a case when the center point is unstable. The parameters are and 
0=1.0 radian. The trajectory passes near a closed separatrix that emerges from and returns to 
the unstable fixed point in the center of the phase space. This description of the motion is 
not exact — in order to have a closed separatrix we need to have a "good" (e.g., analytic) 
integral of motion for the map (12). Such integrals of the motion, in general, do not exist. In 
our case we doubt that the motion is integrable. For large values of the parameter | K\ > 1 
the separatrix for the center point becomes fuzzy and finally the smooth loop-like trajectory 
breaks into a chaotic web around many resonant islands (see Figure 3 for a trajectory with 
the same initial conditions but with /f =-53.3). 




: 



FIG. 2. Unstable behavior of electric field over resonant atomic lattice nodes starting from 

the vicinity of (0,0). Parameter /f =-3.3. 
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FIG. 3. Map (12) trajectory for large absolute values of /f =-53.3 starting from the vicinity of 

(0,0). 

The separatrix in Figure 3 is in reality a very complicated heteroclinic structure, discovered 
by Poincare (see [4] for description of general separatrix motion for symplectic maps). In 
general, there exist outgoing and incoming separatrix curves for any fixed point. These 
curves cross each other with some angle; the crossing point travels away from fixed point 
and back to it in the case when crossing occurs for separatrices of the same resonance. 

The points, traveling from and to the center point are very remarkable — they represent 
a trapped mode of the field that exponentially decays at "+" infinity. In the case of small K 
all the points of the separatrix return back or extremely close to the central fixed point. The 
next section will present an approximation of the map (12) by a differential equation, for 
which the motion turns out to be exactly integrable and the separatrix is a continuous line. 
Appendix A, however, presents a case with random lattice in which separatrix curves don't 
constitute a single line, but there exist crossing points that go from and come to the center 
point, representing a trapped mode with exponential decay of the electromagnetic field at 
both infinities. 



III. APROXIMATE SOLUTION FOR ID STATIONARY DISTRIBUTIONS 

When 0and K are small, the map (12) can be rewritten in the following form: 

4„ + l =q„+ p n 0> 

P , a = ->,,t + P,+^2=- (,4) 
V 1+ 1 i. 1 

Keeping in mind that all changes of p,q are small we can transform the finite difference 
equation (14) into a differential equation, approximating the map: 

ox 

dp , _ Kq (15) 
-JL = -kq + - 
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where k=Q)/ c and / is the distance between lattice nodes. Now we pick the + sign in the 
"source" term and write down the Hamiltonian of the motion (note that coordinate x here 
has meaning of time): 



H 



k(p 2 + q 2 ) K 



(16) 



K 



Only conditions — > k lead to instability and trapped modes: for this condition the central 
point becomes hyperbolic. The equation for the separatrix reads: 



Psep= fkl^ l + ql ' 



(17) 



kl 

where q < 1 . Figure 4 shows the phase space of the Hamiltonian (16) for ft" = 2kl 

K 

One can see two "eyes" with stable motion, a separatrix, corresponding to the trapped 
mode, and an outer region with oscillations around "two-eyed" region. 




FIG. 4. Real (horizontal axis) and imaginary part (vertical axis) of the electric field in ID 
lattice (arbitrary units). The separatrix corresponds to the trapped mode with exponential 
decay of the field at infinity. 

The dependence of the normalized "brightness" (the square of normalized field at the 
center), and the trapped mode size on parameters that can easily be obtained from 
Hamiltonian picture. The normalized "brightness' is roughly the extent of the separatrix in 
the horizontal direction and is: 
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where — > k . The real "brightness" is the normalized brightness multiplied by the 
tiA 

parameter ( ) (see the definition of "normalized field" right before Equation (1 1)). It is 

2ju 

easy to see that the maximum electric field in our trapped modes corresponds to the case 
where the related Rabi frequency equals A. 

The "time" in our case is measured in units of the wavelength and the trapped state size 
is just the time for the field to travel from the maximum coordinate separatrix value to the 
point where the field is reduced by a factor of e. This is, roughly, a half period of oscillation 
around the stable point in the center of the "eye". By finding the stable fixed point, and 
linearizing the motion around it, one can easily find the half-period of the oscillations, which 
is the approximate size S of the trapped state: 



(19) 



where we have — > k again, and the X is the wavelength of the free electromagnetic wave. 

One can see that in the limit K — > kl the size approaches infinity but the "brightness" goes 
to zero. 




IV. TRAPPED MODES FOR 2D AND 3D CASES 



Equation (4) can be generalized to two and three dimensions. Indeed, we can use the 
expression of the spatially averaged source term from (4) in Helmholtz equations. The 
density itself must be small enough to allow the neglect of terms involving the interaction of 
electron dipole moments (e.g., static dipole fields). Definitely, the model is not valid for 
solids - energy shifts of molecules due to electron interactions are too large to ignore. To 
satisfy this requirement, the density N should be less than some parameter K of the order of 
solid density. To derive the equation for 2D and 3D cases, we use here the equation for the 
Hertz vector Z and the equation that relates it to the electric field (see, e.g., [6]): 



V £ Z + k z Z = ±- 



Id 2 

E = WdivZ — T —rZ, 
c dt 

Nk 2 ju 2 . , . . , 

where the new interaction parameter is X = > with JN now being a conventional 

e hA 

atom density. These equations are complicated and rich in various solutions. We would like 
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to reduce the number of cases for analysis by choosing only solutions with div Z — (and, 

consequently, div E — 0). In this case the equations take the form of Equation (4), but for 
vector electric field: 

V z E + k z E = ± 



2* , u2*_^ ZE 



I E I' ' (21) 

div E = 0. 

Kk 

One can see that the ID lattice equations have the same form after substitution — > % . 



The first case we analyze has cylindrical symmetry with the electric field oriented along 
the % axis. The equation for the field E reads: 

d 2 E 1 dE rE 

T + ~ + E z = ~ z (22) 

d(kr) krdkr ' k 2 ^\+\E z \ 2 

where we are interested in only positive % and thus use the plus sign in front of the "source" 
term. Next, we only choose cases for which there is no propagation of waves to 

infinity^ > k 2 c , and we analyze only trapped modes, i.e. the fields exponentially decay at 
infinity. The other criteria for our choice of solutions is that they should be finite at r=0. 
Numerical analysis shows that there exist many (our guess is that there exist infinitely many) 
solutions for trapped modes. They can be enumerated by an integer index that corresponds 
to the number of each solution's zero crossings. Figure 5 shows functions F (left) and F 8 

(right) for the value J /& 2 =1.1 



c As was pointed out in the beginning of the section, the density shouldn't be too large for 
the model to be valid. At the same time, the gas should be dense enough to satisfy a 

necessary trapped mode condition^ > k 2 . Since the parameter depends inversely on 
difference of radiation and transition frequencies, the last condition can always be met. 
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FIG. 5. F„(left) and F g (right) as a function of r. 

One interesting feature of the solutions with large index is that the corresponding trapped 
mode size can be arbitrarily large. It was found numerically for solutions with low index that 
the size and "brightness" of the trapped states are similar to those for the ID case (see 
equations (18) and (19)). Accurate numerical solutions can be easily obtained by simple 
computer programs and are omitted here. Instead, we will construct approximate numerical 
solutions for the trapped modes in the 3D case. 

It is quite probable that the 2D cylindrical mode, if curved into a toroid, produces a 3D 
solution when the aspect ratio of the torus becomes large. To look for similar azimuthally 
symmetric solutions in spherical geometry, we use a particular form of equation (21) for E^, 
independent of <ff. 



d% + 2d_E t + _ L 
dr 2 r dr r 2 sin 



- sin - 



ede 



dG r 2 sin 2 8 h + \E 



(23) 



where 9 is the polar angle. Unfortunately, there is no solution to this equation, independent 
of 6. Therefore, we will look for a solution periodic in#, in the following general form: 



E^r,0) = ^f n (r)P^cos0), 



(24) 



«=i 



where P n (cos 0) are associated Legendre polynomials of index 1 . The choice is obvious - 

the LHS of equation (23) is separable when we substitute Equation (24). Each term (with 
appropriate f„(rj) becomes a solution to LHS of (23) since the polynomials obey the 
following equation: 



1 



sme/- 



sm0d0~"'~ 80 sin 2 6 
The polynomials have the following orthogonality property: 



(25) 
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\pl(t)p;w= 2J ^s nl , (26) 

and the first few of them are: ff (cos 0) = - sin , P\ (cos 0) = -3 sin cos # , 

i 3 2 

P 3 (cos 0) = - — (5 COS 6 — 1) , etc. Substituting (24) into (23), multiplying both sides by an 

associated Legendre polynomial with index /, and integrating the resulting equation using 
condition (26) over 6, yields: 



^ + H +/i( ^-W ) = |g)yj W^cos^sinfl 

\l + (Zfn(r)P:(cos0)) 2 

V n=l 

As in the case of cylindrical solutions, we found many approximate solutions for some small 
harmonic numbers. One of these harmonic numerical solutions is shown in Figure 6. The 
figure shows the cross section of the trapped mode when y=0. One can see that it consists 
of two torroids (rings) one atop of another. Of course, finding numerical solutions is no 
proof that the procedure converges as the number of harmonics increases, but this approach 
works in most cases and we assume this is the case for our equations. 

We guess that more spherical forms of solutions can be obtained when we try to find 
solutions for and Eg, rather than E^. Unfortunately, in this case the equations become 
more involved and the analysis is left for future developments of the topic. 




FIG 6. Electric field squared in 3D trapped azimuthally symmetric mode as a 
function of X, Z coordinates. Shown is the cross section of the mode with Y 
coordinate equal to zero. 

In addition to the stationary distribution, one may look for traveling distributions. It is 
straightforward to modify equations (4) into traveling wave equations. Once again, we leave 
this subject to future exploration of the topic. 
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V. MODEL LIMITATIONS AND PRACTICAL EXAMPLES 



Throughout this paper we have neglected the spontaneous emission of radiation, or, in 
general, the quantization of the electromagnetic field. As was shown in some papers on Rabi 
oscillations (see, e.g., [7]), the semiclassical theory leads to almost the same qualitative 
predictions as the exact theory with field quantization, including the correct rates of 
spontaneous emission. The trapped mode solutions will be the same in the exact model in 
the limit of high stored energies in the trapped modes. In other words, the spontaneous 
radiation will be also trapped as well as the stimulated emission - there is no energy escape 
from the trapped modes even when the fields are quantized. When atomic motion is slow, 
and the ionization and other energy consuming effects are small, the fields can be trapped 
for very long times. 

One more important physical property of the trapped modes should be discussed before we 
give estimations for available fields. Obviously, trapped fields exert pressure on the resonant 
atoms similar to the RF cavity case. To counteract the field pressure, we have to consider 
atomic restraining mechanisms in the medium. For example, the ID atomic planes can be 
squeezed in between transparent solid slabs, or atoms can be placed in some crystal lattice in 
analogy with photonic crystals, forming some (not necessarily periodic) lattice. If the atoms 
are in a gas and the pressure from the fields is comparable with the gas pressure, in addition 
to given equation we must also solve the equation for pressure (or atomic density) that 
results from magnetic fields (since the magnetic field in our solutions is always in phase with 
currents). For steady state solutions with zero average velocities the equation reads: 

VP = f, (27) 



where P is the pressure and/is the magnetic force, acting on the atomic currents. Obviously, 
in strong standing electromagnetic waves the atoms in the gas will form a diffraction grating 
for resonant atoms that will lead to even stronger trapping of waves. The total self consistent 
equation for the atomic density and the resonant fields will be treated in future explorations 
of this topic. 

In principle, this electromagnetic field trap is similar to superconducting RF cavities. If 
we form (even temporarily) structures of standing or slowly moving atoms similar to RF 
cavities, the trapped modes can be used to provide acceleration. We now estimate the 
achievable fields in the trapped modes. The rotating wave approximation (RWA), used 
above, fails (see [8]) when the Rabi frequency becomes comparable with the transition 
frequency between two levels in our approximation, and, considering multi photon effects, 
the Rabi frequency should be less than the frequency difference between our chosen level 
and all closest levels. We take as an example the hydrogen atom and its levels with n=l and 

2 8 ea 

n=2. The transition frequency is GJ^ ~l.55ld 6 Hz and jU 12 = -\d 3 r u[{r)ezu 2 (r) = " 

J 3 V2 

for the transition between the 1 st and 2 nd states, where the upper level has L=l, and m=0 

quantum numbers (the reference axis here is taken to be the direction of electric field), and 

a is the Bohr radius. Assuming the Rabi frequency is only 1% of the transition frequency, 

the distance between levels n—2 and n—3 is roughly fourteen times more than the Rabi 

frequency. Therefore, our approximation is reasonable for such values of the Rabi 
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frequency. Equating the Rabi frequency Q. = jLlE Q ITi to one hundredth of the transition 
frequency, we get a maximum amplitude (see the elaboration on maximum electric field after 
Equation (1 8)) of electric field E ~2.5 GeV/ m, that is almost 2 orders of magnitude larger 
than that of the achievable electric fields in superconducting cavities. 

Unfortunately, the density of the Hydrogen gas must be very large. The trapped mode 

Nk 2 U 2 2 
parameter^ = must be larger than k . After equating these and substituting the 

e Q hA 

2 ea n 

above jU = jU l2 = - ,— into the expression for the parameter^ we obtain a minimum 

3 V2 

hydrogen density of N-3.3-10 21 cm 3 . This density is achievable at 100 atmospheres pressure. 
For some metals, the dipole transition coefficients can be an order of magnitude larger than 
that of hydrogen, therefore these fields can be achieved for pressures comparable with that 
of the atmosphere. For hydrogen at atmospheric pressure we obtain achievable fields of 25 
MeV/ m, which is also not small. But in any case, the total pressure of the fields can't be 
larger than the atmospheric pressure. For trapped modes with the radius R=20 cm the total 
energy of the state can be as large as P-47tR 2 =10kJ. 

We have neglected all energy consuming mechanisms, such as gas heating and gas 
ionization. However, for ultraviolet light, the air breakdown for wavelength of 0.1 (Im 
requires a field of 10 9 -10 1() V/m ([9]). Therefore, the estimated fields can exist in the air 
without ionizing the gas. This model is one possibility to explain ball lightning - it may be a 
trapped electromagnetic wave whose frequency is close to the transition frequency between 
some atomic levels of the air atoms. In addition to trapping, some mechanism should exist 
to transfer energy from, for example, discharge in the air, to the trapped wave. These 
questions are beyond the scope of the present paper. 



VI. CONCLUSION 

This paper demonstrates that there exist a variety of trapped solutions for 
electromagnetic fields in a medium where resonant atoms are present. The total reflection of 
electromagnetic waves by the medium is a necessary condition for their appearance. In 
addition to total reflection, the nonlinearity of atom-wave interaction is responsible for the 
phenomenon. Namely, for large intensities of the electromagnetic field, the "source term" 
(see RHS of Equations (15) and (20)) saturates and the medium becomes transparent to the 
wave. The region of "transparency" forms the core of the trapped mode. At the edges the 
field vanishes and the medium reflects the wave completely, thus confining the localized 
state with high electromagnetic energy density. 

The trapped modes, in principle, may have zero loss and accumulate enormous field 
energy. It may be possible to use trapped modes in the same way as RF cavities or wave 
guides are used, but with larger achievable fields and with potentially zero energy loss. In 
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addition, the trapped solutions for electromagnetic fields may be responsible for some 
atmospheric phenomena, such as "ball lightning". 
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APPENDIX A. TRAPPED MODES IN RANDOM ID LATTICE 

Localized states also exist in random lattices. In addition to mode trapping, which is related 
to nonlinearity of field-atom interaction in the semiclassical approach, the randomness adds 
one more effect, known as "localization" in solid state physics [2]. We do not present a 
complete review the phenomenon here. However, we do point out one mathematical 
theorem by Furstenberg [5], which states that for any product of a large number of random 
noncommuting ID (2x2) matrices with unit determinant, one of the eigenvalues of the 
resulting matrix will grow exponentially with the number of multiplication terms (the other 
eigenvalue is, of course, the inverse of the growing one). In application to our problem, a ID 
random lattice of atoms will always yield exponentially growing (decaying) fields at infinity. 
Following Poincare (see, e.g. [4]), we define the sequence of points that decay exponentially 
at — oo as the H" curve, and the sequence of points, decaying at +°°, as H + , respectively. 
These curves are not trajectories, they are collections of trajectories leaving and returning to 
the center of phase space at infinities. In our case we have two H + and H curves, each curve 
having its symmetric counterpart with respect to 180 degrees rotation of the phase space 
around the center. They are related to each other only by the fact that both have the same 
asymptotic growth/ decay rate. The H" curve, in general, grows to substantial values 
(probably to infinity) at +°°, as seen in numerical examples, but this growth is slow because 
the kick from each atomic node becomes constant for large values of the field (in our units, 
large as compared to unity). The H + curve behaves similarly, at -oo. These curves, in general, 
are not the same for periodic systems, with the "small" exception of integrable systems. 
When the phase space, where these curves evolve, is confined or the growth is slow, as in 
our case, the curves cross each other with some nonzero angle (probably, there exist some 
exceptions, but we haven't seen them in our numerical cases). The crossing point is our 
trapped mode — it grows from zero values at — °° and decays exponentially at +°°, because it 
belongs simultaneously to both "+" and "-"curves. The crossing of curves is confirmed by 
all numerical cases performed by the author (that, of course, does not constitute a 
mathematical proof). Here, we present one example of a random lattice with all the 
described curve and trajectory species, selected randomly, from infinitely many random 
cases. 

For this case, we choose K=-0.8, 100 random nodes and random phase advances from to 
1 radian (chosen as a fractional part of the node number after it is divided by 71) between 
nodes. The approximate H" line was built in the following way: we selected a very small 
initial value of the forward field (see u f in equation (11)), evolved its growth through 10 
nodes and connected first and last point by a line. Thus, we approximated the separatrix, 
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coming emerging from zero for — °°. We approximated this line by 1000 points placed 
uniformly over the line extent and tracked it over 50 nodes. The solid line in Figure 4 shows 
the resulting line in the "phase space" of the map, where the x coordinate stands for the 
imaginary part, andj - for the real part of the field. This line is an approximate piece of the 
H" line. The same was done for the inverse map, starting from +°°, that corresponds to the 
100" node, back to the center of the lattice (node 50). The dashed line in Figure 7 represents 
part of the H + curve, obtained in this approximation. One can see that the curves do cross 
each other. 



Fig. 7. Separatrix crossing. 

The crossing point tracked from the 1 st to the 100 th node is shown as a trajectory in Figure 
8 (dashed line). It starts at -2.6*10" 44.0*10 forward field value at node 100 and ends near 
zero at the 1 st node (the actual values obtained are at 10"" level, due to the limited accuracy of 
the numerical simulations and the strong influence of errors in the case with exponential 
instability). The solid line shows the trajectory of the H" family that passes very close (within 
0.01 at node 50) to the crossing point. Due to a slight error in initial conditions and the 
exponential instability, this trajectory eventually diverges from that of the trapped mode. 



FIG. 8. Trapped mode shown by the dashed line. The horizontal coordinate corresponds to 
the node number, and the vertical coordinate shows the modulus of the forward wave. The 
solid line shows the field behavior for initial conditions 10% different from those of the 
dashed line. 
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